Astronomy & Astrophysics manuscript no. mhdfil 


©ESO 2012 


July 17, 2012 





Profiles of interstellar cloud filaments. 
Observational effects in syntlietic sub-millimetre observations. 

M. Juvela^ , J. Malinen^ , T. Lunttila^ 



O 

(N 

O: 
o 



Department of Physics, P.O.Box 64, FI-00014, University of Helsinki, Finland, mika.juvela@helsinki.fi 
Preprint online version: July 17, 2012 

ABSTRACT 

Context. Sub-millimetre observations suggest that the filaments of interstellar clouds have rather uniform widths and can be described 
with the so-called Plummer profiles. The shapes of the filament profiles are believed to be closely linked to the formation mechanism 
and physical state of the filaments. 

Aims. Before drawing conclusions on the observed column density profiles, we must evaluate the observational uncertainties. We 
want to estimate the bias that could result from radiative transfer effects and from the variations of submm dust emissivity that are 
expected to take place in dense clouds. 

Methods. We use cloud models obtained with magnetohydrodynamic simulations and carry out radiative transfer calculations to 
produce maps of sub-millimetre surface brightness. Column densities are estimated based on the synthetic observations. For selected 
filaments, the estimated profiles are compared to those derived from the original column density. Possible effects from spatial varia- 
tions of dust properties are examined. 

Results. With instrumental noise typical of the Herschel observations, the parameters derived for nearby clouds are correct to within 
a few percent. The radiative transfer effects have only a minor effect on the results. If the signal-to-noise ratio is degraded by a factor 
of four, the errors become significant and for half of the examined filaments the values remain practically unconstrained. The errors 
also increase in proportion to the cloud distance. Assuming the resolution of Herschel instruments, the model filaments are barely 
resolved at a distance of ~400pc and the errors in the parameters of the Plummer function are several tens of per cent. 
Conclusions. The Plummer parameters, in particular the power-law exponent p, are sensitive to noise but can be determined with 
good accuracy using Herschel data. However, one must be cautious about possible line-of-sight confusion. In our model cloud, a 
large fraction of the filaments identified from column density maps are not continuous structures in three dimensions. 

Key words. ISM: clouds - Infrared: ISM - Radiative transfer - Submillimeter: ISM 
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1. Introduction 

For a long time the filaments have been recognised as a com- 
mon feature of interstellar clouds that is also closely related 
to the star formation process (e.g . . lElm egreen & El megreenI 
Il979t ISchneider & Elmegreenl l 19791 'Bai Tv et alJ ll987). Recent 
Herschel studies and ground based surveys have revealed 
the ubiquity of filamentary structures within dense interstellar 
clouds, and pre-stellar cores and protostars along dense, grav- 
itationally unstable filaments (e.g., [M en'shchikov et al J l201(]t 
lAndre et alJl2010HArzoumanian et al.ll201 It lJuvela et al. l2012h . 
Therefore, the structure and formation of the filaments are 
important parts of star formation studies. The filaments are 
naturally form ed as a result of interstellar turbulence (e.g., 
iPadoan & No rdlund 201 1), but significant contribution must ex- 
ist from immediate triggering by supernova explosions and 
the radiation and stellar winds from massive stars (e.g., 
iPeretto et al]l2012l) . The filaments are expected to undergo grav- 
itational fragmentation that will lead to star formation (e.g., 
[inutsuka & Miyama 1997; Myers 2009; Heitsch et al. 201 1). 

Many recent studies have addressed the question of the typ- 
ical shape of the filament cross sections. The observed col- 
umn density profiles are fitted with the so-called Plummer func- 
tion, partly because of the theoretical predictions associated with 
this functional form. In particular, an infinite isothermal cylin- 
der in hydrostatic equilibrium is exp ected to follow this profile 
(IStod6Lkiewicj[T963l : IOstrikeilll964t) while, for example, mag- 



netic support, deviations from the isothermal conditions, and ini- 
tiated collapse would all cause char acteristic effects on the pa - 
rameters of the profile function (see lArzoumanian et aL I I2OII . 
If the parameters of the profile function can be determined with 
high accuracy, this will provide important information on the for- 
mation and the current state of the filaments. 

The central question is, can we accurately measure the struc- 
ture and mass of the filaments. Our knowledge is mostly based 
on molecular lines and dust emission. The kinematic information 
provided by line observations is crucial but the use of molecules 
for a structural analysis is complicated by the abundance and ex- 
citation variations that are difficult to measure and to model. The 
situation is exacerbated in the dense environments where many 
molecules are severely depleted as they begin to freeze onto dust 
grains. 

Dust emission is considered a reliable probe of dense clouds, 
especially at sub-millimetre wavelengths where the dependence 
on the dust temperature is reduced. However, the line-of-sight 
temperature variations can still play a significant role and, in the 
case of dense cores, can rend er the column densi ty estimates un- 
certain by a factor of several dMalinen et alJ20fl ) . Furthermore, 
dust is expected to undergo evolution in dense and cold envi- 
ronments. This may include accretion of ice mantles, agglom- 
eration of dust particles (leading to the disappearance of small 
grains), and temper ature-dependent changes in t h e dust opti- 
cal properties (e.g., lOssenkopf & HenningI Il994t lOrmel et aP 
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Table 1. Main parameters of the model clouds. 



Model 


<Av> 
(mag) 


< il(H2) > 

(cm-3) 


Dust model 


A 


1.0 


111.0 


standard MWD' 


B 


2.0 


222.0 


standard MWD 


C 


4.0 


450.0 


standard MWD + dust with higher 








sub-mm opacity 



' Milky Way dust with Ry = 3.1. 



1201 itlMenv et al.l2007l:lKohler et al.l20T J!). The analysis of sub- 
millimetre observations is particularly affected by changes in the 
dust spectral index /3 and by the possible increase of the dust 
opacity k. The value of k may vary by a factor of several be- 
tween diffuse and dense regions. This is important for the study 
of cloud structure, because the derived column densities are di- 
rectly proportional to 1/k. These effects could be expected to be 
less severe in filaments than in the cloud cores because of the 
lower density of the filaments and the consequently longer time 
scale of grain growth. However, significant changes in dust sub- 
millimetre opacity have been seen also towards filaments. For 
example, a factor of -3.4 i ncrease was reporte d for a filament 
in Taurus with Ay < 10™ (IStepnik et al.ll2003h . In addition to 
the uncertainty of k, a wrong assumption of /? will bias the mass 
estimates throiigh its effect on the colour temperature estimat es 
(iMalinen et al.ll201 ll: LTuvela & Ysardll2012l: lYsard et al.ll2012h . 

In this paper we investigate the possible differences between 
the real density structure of the filaments and the one deduced 
from typical sub-millimetre observations of dust emission. The 
cloud models are based on magnetohydrodynamic simulations 
with distinct filamentary structures. The emphasis is not on the 
nature of the filaments themselves but on the uncertainties of de- 
riving their properties from observations. We carry out radiative 
transfer modelling to produce synthetic surface brightness maps 
and estimate the profile functions of selected filaments. The re- 
sults are compared to the real column density structure read di- 
rectly from the model. 

The structure of the paper is the following. In Sect. |2] 
we present the cloud models and the calculations leading to 
synthetic surface brightness maps. In Sect. [3] we describe the 
selected filaments and the procedures used to fit them with 
Plummer profiles. The main results are shown in Sect. |4] where 
we list the results for the filament profile fits, both using the true 
column density and the column density estimated from the syn- 
thetic observations. We discuss the results in Sect.|5]before pre- 
senting our final conclusions in Sect. |6] 



2. Simulated clouds and observations 

The modelling is based on a magnetohydrodynamic (MHD) sim- 
ulation carried out on a regular grid of 1000^ cells. The rms 
sonic Mach number is ~9 and the rms Alfvenic Mach number 
~2.8. The data are the same that were used in iMaUnen et alJ 
(1201 Ih and we refer to that paper for details. The MHD calcu- 
lations assume an isothermal equation of state. The simulation 
was run first without gravity and the stirring continued for sev- 
eral shock crossing times to develop turbulence well separated 
from the initial conditions. The simulation was continued with 
gravity switched on and the data used here corresponds to the 
situation during early co llapse (Mod el I in Malinen et al. 201 1). 
We use the unigrid run dPadoan & N ordlund 2 01 11) rather than 
the Adaptive Mesh Refinement (AMR) MHD data because we 



need high spatial resolution not only for the cores but also for 
the lower density filaments. 

The cloud size was scaled so that the mean visual extinction 
Ay through the cloud is l.O" or 2™. In the following, these are 
called models A and B. Assuming a linear cloud size of 6pc, 
these coiTespond to mean hydrogen densities of 1 1 1 cm"^ and 
222 cm"-', respectively. 

For the radiative transfer calculations we built an adaptive 
grid from the 1000^ full resolution model. The adaptive grid is 
constructed by combining neighbouring cells that are optically 
thin, and have similar densities. The maximum allowed V band 
optical thickness of a cell is t^^^ - 0.04 in model A. The max- 
imum density contrast was set to Pmax/Pmin = 30. The resulting 
adaptive model has approximately 160 x 10'' cells and consists 
of 839 grids. The same grid was used for model B. The grid- 
ding ensures that also moderately dense region, the filaments in 
particular, are sampled at the best resolution. 

We use the radiati ve transfer program described in 
Lunttila & JuveH (1201 2|) (see also lJuvela & PadoanI 120031; 
Juvela 2005). The cloud is illuminated externally by the isotropic 
interstellar radiation field (Mathis et al. 1983). The dust is as- 
sumed to have the properties of the stan dard Milky W ay diffuse 
interstellar dust (Rv-3.l) as described in lDrain^ ( 12003 ). Because 
of the low dust temperatures and the modest optical thickness of 
the model clouds, the absorption of the thermal dust emission 
does not contribute significantly to the dust heating. Therefore 
the dust temperature was calculated in a single step without iter- 
ation. 

In the third model C the average visual extinction is set to 
A V ~ 4™. Assuming a linear cloud size of 6 pc, the mean density 
is n(H)=450cm"^. In model C we use two dust components. The 
second one coiTesponds to a modification where the dust sub- 
millimetre opacity is increased. This is implemented by multi- 
plying the original opacities with a factor (/l/20/im)°-^ at wave- 
lengths longer than 20 //m. The resulting dust opacity spectral 
index is thus /? ~ 1 .5 in the sub-millimetre and, as a result, the 
sub-millimetre opacity Ky is higher For example at 250 fim the 
Ky has increased by a factor of 3.5. The relative abundance of the 
two dust components is varied according to the local density so 
that the abund ances are equal at a density of n(H) ~ 5 x lO'* cm"^ 
(see Fig. 5 in iMalinen et al.ll201 lb . At low densities the relative 
abundance of the modified dust component goes to zero. By the 
time the density reaches 10^ cm"-', practically all the dust is in 
the modified form. The calculations for model C were canied 
out at a resolution of 500^* cells where the grid was obtained by 
direct downsampling of the orig inal 1000^ data. 

The radiative transf er calculations ( lJuvela & PadoanI 120031 : 
iLunttila & Juvelall2012l) solve the 3D temperature structure of 
the model clouds and result in synthetic surface brightness maps 
of 1000x1000 (models A and B) or 500x500 pixels (model C). 
Surface brightness maps were calculated at 160, 250, 350, and 
500 jum that coiTespond to the wavelengths of Herschel observa- 
tions. The full-resolution images were convolved with Gaussian 
beams with FWHM corresponding to 3, 6, or 12 pixels (in model 
C the beam widths in pixels are half of these values). Before the 
convolution we add to the maps white noise that corresponds 
to 3.7, 1.20, 0.85, 0.35MJy/sr per original Herschel beam size 
of 12, 18, 25, and 37", respectively. The values were obtained 
using the Herschel Interactive Processing Environment HIPl^l 
and they correspond to the instrumental noise levels of typi - 
cal Herschel mapping observations (as in Mali nen et al.ll201 ll) . 
Most actual Herschel observations are confusion Umited, partic- 
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Fig. 1. Ay maps of the selected twelve filaments in model A. The dashed lines indicate the filament segments included in the 
analysis. The Ay values are obtained directly from the model cloud where the 1000x1000 pixel image has been convolved with a 
beam with FWHM equal to three pixels. 



ularly at wavelengths A > 250 fim. In our simulations the con- 
fusion results only from structures inside the modelled volume. 
In more realistic simulations, especially for low surface bright- 
ness values, one might also have to consider the confusion by 
extragalactic point sources. The final common resolution of all 
the maps (i.e., the size of the beam used in the convolution) is 
40". When the data are convolved with a beam equal to three 
pixels, the pixel size is thus 13.3" and the map size is 3.7 de- 
grees. With the assumed cloud size of 6pc, the FWHM beam 
sizes of 3, 6, and 12 pixels translate to the cloud distances of 
93, 186, and 371 pc (in the following these are called the near, 
medium, and far distances). For the near distance, the filaments 
are clearly larger than the beam. For the medium distance, the 
filaments are still slightly larger than the beam, but in the case 
of the far distance, the observed filament width is dominated by 
the convolution. 



3. Methods of analysis 

The synthetic surface brightness observations were converted to 
maps of column density. The dust colour temperature was cal- 
culated pixel by pixel by fitting the observed 160, 250, 350, and 
500 yum intensities with a modified black body curve. The data 
points were weighted according to the noise levels of the simu- 
lation. In the fits the opacity spectral index was fixed to a con- 
stant value of j6 = 1.8 or/3 = 2.0. These values differ from the 
actual value in the dust model (MiUcy Way dust with Ry - 3.1) 
where the spectral index between the wavelengths of 160 fim and 
500 yum is 2.09. 



With the dust colour temperature known, the column density 
is calculated as 



N' 



Obs 



ByiTc)Ky 



(1) 



where 5 y is the intensity of the fitted modified black body curve. 
By is the Planck function, and Ky the dust opacity. The Ky value 
was taken directly from the dust model. Therefore, the derived 
column density A^°''^ diff'ers from the true column density A^t™'= 
only because of noise and because of systematic errors in the 
colour temperature Tc, the latter being caused by line-of-sight 
temperature variations. 

We selected four filaments from each of the three column 
density images calculated towards three orthogonal directions 
(see Fig.[TJ. The filament segments were selected by eye and the 
path of each filament was set to follow the peak of the filament 
as seen in the A^^™ maps obtained from the cloud model. 

The pr ocedure used in an alysing the filaments is basically the 
same as in lJuvela et al.l (12012 ). The column density profiles per- 
pendicular to the filament path were extracted at 20" intervals. 
The individual profiles were aligned so that their peaks match 
before calculating the median profile of the whole filament. Thus 
the location of the filament was allowed to differ from the orig- 
inally selected path. A fixed range of [-0.3 pc, -1-0.3 pc] around 
the peak of the column density profile was included in the sub- 
sequent analysis. The filament FWHM was measured by fitting 
the profile with a Gaussian function plus a constant background. 
The filament FWHM values were determined for each cross sec- 
tion separately and for the final combined median profile. The 
median profile of each filament was also fitted with a model 
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Fig. 2. The profile of the filament number 1 as seen in the models 
A, B, and C (frames a, b, and c, respectively). The solid thick 
line shows the actual column density profile in the model cloud. 
The other curves show the profiles derived from the synthetic 
observations when the cloud is at a near distance of 93 pc (thin 
solid line) or 371 pc (dotted line). 



consisting of a linear background plus a Plummer-Uke column 
density profile 



N(r) = A„ 



Pc^flat 



(2) 



'[l+(r//?flat)2](''-l'/2 

/see lArzoumanian et all 1201 it iNutter et al.1 120081: ljuvela et alJ 



|2012|). The factor A 



rr(i + MV /^ cosn 



is calculated using the formula Ap - 
du, where we assume an inclination an- 
gle of ; = (lArzoum anian et al. 2011). In addition to the two 
parameters describing the background, the free parameters of 
the Plummer function are the central density of the filament p^, 
the size of the inner flat part TJflat, and the power-law exponent 
p. The finite resolution of the observations was taken into ac- 
count by fitting the observed profile not directly with the theo- 
retical Plummer profile but to a Plummer profile convolved with 
the beam used in the simulated observations. The convolution 
weights take into account the fact that the actual convolution 
takes place over a 2D map. This allows us to derive deconvolved 
parameters for even unresolved filaments. 

In addition to the fitted parameters, we also look at some 
derived quantities. In particular, the mass per unit length of fil- 
ament, Miine, can be calculated by integrating over the fitted 
Plummer profile. 

4. Results 

The twelve filaments were analysed using the true column den- 
sity A^T'™ and the column densities N'^^'^ estimated from syn- 
thetic observations. The analysis was repeated with an actual 



resolution equal to 3, 6, and 12 pixels of the original 1000x1000 
pixel maps, these corresponding to the near, medium, and far 
cloud distances. Because of the lower resolution of the model C, 
these correspond there to convolution with a beam with FWHM 
equal to 1.5, 3, and 6 pixels, respectively. 

Figure |2] shows examples of the column density profiles, 
comparing the true column density profiles to those estimated 
from the synthetic observations. These already indicate that in 
models A and B the column densities are underestimated by 
~10% at the centre of the filament, because the estimated colour 
temperature Tc overestimates the mass-averaged dust tempera- 
ture. The fact that the analysis is carried out with a slightly too 
low value of the spectral index, /3 - 2.0, instead of 2.09, con- 
tributes in the same direction but is a smaller factor In model 
C, the sub-millimetre dust opacity was increased by a factor of a 
few in the densest parts of the cloud. This is directly reflected in 
the column density estimates that are more than twice the correct 
value at the centre of the filaments. The figure also shows that, 
apart from the filament itself, the profiles are affected by other 
emission components, i.e., by background confusion. 

Each filament was fitted with a Plummer function. Fig. [3] 
shows the results for model A in the case where the observed 
column densities were estimated with /3 = 2.0 and the cloud 
distance was set to d - 93 pc, which corresponds to convolu- 
tion with a beam with the FWHM equal to three pixels. The 
figure demonstrates that the profiles can be well fitted using the 
Plumm er function. In Fig. [3] the parameter values show a large 
scatter. iMalinen et al.l (12012) noted that the fitting of Plummer 
profiles is sometimes sensitive to the initial parameter values of 
the least squares minimisation. Visual inspection shows that all 
the fits in Fig.[3]are good. However, in addition to intrinsic vari- 
ations between the filaments, some of the scatter is caused by 
the instrumental noise and especially by the confusion caused 
by other emission along the line-of-sight. 

Figure ID shows a comparison of the estimated values of p^, 
p, and Miine for all three beam sizes (3, 6, and 12 pixels) for 
model A. On the x axis are the parameters obtained using the 
true column densities, on the y axis are the values derived from 
the analysis of synthetic observations with and without obser- 
vational noise and with /3 equal to 1.8 or 2.0. The effect of the 
A/3 = 0.2 is relatively small and is seen mainly as a systematic 
effect in the density pc and mass Mum values where it amounts 
to an underestimation of 20-30%. Without noise the parameter 
values obtained using observed and true column density are sim- 
ilar to within a few percent. The effect of noise is often small but 
can also lead to tens of percent difference in the parameter val- 
ues. The loss of resolution (i.e., convolution with a larger beam) 
has a clear impact when noise is included and leads to larger 
scatter especially in the estimates of the central density, which 
are underestimated in most cases, even up to ~2 orders of mag- 
nitude. In the fitting routine, the maximum allowed value for the 
parameter p is 10. The largest bias in p values occurs when the 
value derived from either the true or observed column density 
(but not both) is ~10. Compared to the other parameters (espe- 
cially the central density pc), the mass MLine is less affected by 
the noise and the resolution. 

The corresponding results for the higher column density 
model B are presented in Fig.|5] Also here the parameter values 
obtained using observed and true column density without instru- 
mental noise are in most cases similar to within a few percent. 

The scatter of the central density pc values (from A^°''^ in- 
cluding noise vs. those from N^""^) appears to be lower com- 
pared to model A. Using the observations with beam sizes equal 
to 3 and 6 pixels (corresponding to the near and medium cloud 
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Fig. 3. Fits of the Plummer function to profiles extracted from the observed column density maps (Ji-l.Q) of Model A cloud at 
the near distance. The solid black curves show the estimated median profiles of the filaments and the dashed cyan lines the fitted 
functions. The parameters pc (central number density of H2 in units of cm"^), p, and Miine (in units of Mq/pc) are shown in the 
frames as well as the rms value of the fit (in units of 10^' cm"^) and the FWHM values of the filaments obtained from fits with 
Gaussian functions. 



distances) the density and mass are mostly either correct or 
somewhat underestimated. However, when low resolution with 
a beam size of 12 pixels is used, the scatter increases and the 
parameters can be also overestimated. For the Mune values the 
error is up to a factor of ~2. One should note that this is the error 
between the masses derived from either the true or the estimated 
column density. We will return later to the comparison with the 
actual mass of the filaments that can be obtained directly from 
the model cloud. 

In model C, the abundance of modified dust, i.e., the compo- 
nent with higher sub-millimetre opacity, increases with density. 
Figure |6] shows the relative abundance of the modified dust as 
projected values towards three orthogonal directions. The loca- 
tions of the selected filaments are marked in the figure. The fig- 
ure shows that while modified dust is the dominant component 
towards the centre of many filaments, its abundance remains be- 
low 20% in most of the area included in the Plummer fits, i.e., 
at distances of a few arcmin from the filament centre. Thus, the 
spatial variation of dust properties should be able to bias both 
the A^obs values and the observed shape of the profiles and thus 
affect the parameters of the Plummer fits. 

Figure [7] presents some examples of the column density 
profiles of the model C filaments, comparing 1^"°^ and N"^™'^ . 
Because the column densities are estimated assuming a constant 



dust opacity, Ky, the column densities are overestimated when- 
ever the dust opacity has been increased. In the figure the N"^™'^ 
values have been scaled by a factor of 1.5. This provides rough 
agreement with N^'°^ in the outer parts of the filaments, while at 
the centre the synthetic observations lead to still higher column 
densities by a factor of ~2. For example, the effect on the ob- 
served filament FWHM depends on whether the emission at the 
half value point is produced mainly by the normal dust or the 
modified dust. In the latter case, the observed filament FWHM 
would not be affected even though the column density can be 
strongly overestimated. 

The values of the Plummer parameters obtained with model 
C are shown in Fig. [8] The optical depth of the cloud is about 
twice as high as in the case of model B. Therefore, in addition to 
the spatial variation of dust properties, also the radiative transfer 
effects should have a stronger effect on the results. The observa- 
tions are seen to overestimate the filament mass up to a factor 
of several. Thus, the increase of the average dominates over 
the opposite effect resulting from the dispersion of dust tempera- 
tures along the line of sight. The volume density estimates can be 
~0.5-l orders of magnitude too high and are more biased when 
the resolution of the observations is low. The scatter is larger 
than in the previous models, especially when the data are con- 
volved with the largest 12 pixel beam. The smaller spectral index 
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Fig. 4. Comparison of the Plummer parameters derived from true column density and the values estimated from the synthetic maps 
of model A. The rows correspond to parameters p^, p, and Mune- The columns correspond to beam sizes equal to 3, 6, and 12 pixels, 
i.e., to the near, medium and far cloud distances. On the x-axis are always the parameter values from A^t™'= and on the y-axis are the 
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Fig. 7. Profiles of the filaments in model C (at the near cloud 
distance). The thick lines correspond to the true column density 
^Trae q£ jjjg jjjodel cloud multiplied by a factor of 1.5. The thin 
lines denote the column density derived from synthetic observa- 
tions. 



/3- 1 .8 still leads to smaller MLine values but the scatter of the val- 
ues suggests that the fits are less well constrained. The p values 
can be divided to two groups: those close to p ~IQ and those 
with values p ~2-4. At low p the values obtained from synthetic 
observations do not show any systematic bias. Contrary to mod- 
els A and B, in model C there are also cases without instrumental 
noise where the p value derived from either the observed column 
density or the true column density (but not both) is close to 10. 

The FWHM values of the filaments were calculated with 
Gaussian fits and Fig. |9] compares the values in the three mod- 
els. The values derived from A^t™'= should be identical in all the 
three models. However, because of the coarser initial resolution 
of model C, the values of that model differ slightly. In models 
A and B the values obtained from synthetic observations are al- 
ways within a couple of per cent of the A^t™'= results. In model 
B the filament optical depths are above Ay - 10™ and one could 
expect the central column density to be underestimated relative 
to the lower density outer regions. However, the effect on the 
beam convolution is still almost negligible. More dispersion is 
seen in model C. Because the enhanced emissivity increases the 
signal at the centre of the filaments, the estimated FWHM val- 
ues of the filaments would be expected to decrease. This is also 
the observed trend, the analysis of A^'-''"* leading to ~ 10 - 20% 
lower values than A^^™*^. The effect is not quite constant, prob- 
ably because of the additional effect of instrumental noise and 
the confusion by dense background structures. The largest dif- 
ferences are ~40% (filaments number 1 and 12). 
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Fig. 5. As Fig. |4]but for the model B that has twice as high column density. 
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Fig. 6. The fractional abundance of modified dust in the model C. The three frames correspond to three orthogonal directions of 
observations. The white circles indicate the ends of the filament segments included in the analysis. 



We compare the p values of the twelve filaments in the three 
models in Fig. [TOl using the medium cloud distance of 186pc 
(beam size equal to 6 pixels) and /? - 2.0. When the parameter 
values are derived from the true column density, A^^™*^, (frame 
c), all the models give similar results regardless of the value. In 
four out of twelve cases the p value is above 8, indicating that the 
parameters are not well constrained, i.e., the minimum value 
is reached with parameter values that are far from those obtained 
in the absence of noise. Most of the other twelve filaments have 
p values in the range 2^. When the parameter values are de- 
rived from the observed A^*-""* without instrumental noise (frame 
a), the results are mostly similar (to within Ap ~ 0.5) to the pa- 



rameter values derived from the true . Only in two cases 
the p value in model C is biased by a factor of ~2, notably to op- 
posite directions. However, when using the observed A^*^''^ with 
noise (frame b), the situation gets more complicated. Now only 
half of the filaments have similar p values (to within Ap ~ 0.5) 
as in the case with the true N^"^'^. This applies to all three mod- 
els. The bias can now also be seen in all models, in particular 
in model A. These results indicate that the effect of instrumen- 
tal noise on the derived p parameter is more prominent than the 
radiative transfer eff'ects. 

The results of Figs.lH |5j and |8] are summarised in Fig.fTTI 
The figure shows for the three model clouds and the three cloud 
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Fig. 8. As Fig.|4]but for the model C containing two dust components with a spatially varying abundance ratio. The beam size (3, 
6, or 12 pixels corresponding to the near, medium, or far cloud distance, respectively) is indicated on top of each column of frames. 



distances the distributions of the relative parameters errors using 
boxplotfl The errors are calculated as the difference between 
the values derived from A^'-"'* (with fi-l.Q) and N"^™'^, normalised 
with the latter parameter estimate. The figure includes the error 
distributions for the logarithm of the central density pc, parame- 
ter p, and the mass Mune- 



5. Discussion 

The estimated central density of the filament pc, the power-law 
exponent p, and the mass per length unit MLine were examined in 
Sect.m We have shown that the filaments of our MHD model can 
be well fitted with the Plummer function. However, the fitting is 
sometimes sensitive to the initial values of the fitting parame- 
ters (not all initial values leading to a good fit) and the param- 
eters extracted from apparently good fits can have a large scat- 
ter, especially in the case of filaments that are not well resolved. 
This indicates that the parameters are sometimes sensitive to the 
amount of instrumental noise added and thus very likely also 
to the confusion by background structures. All the fits used in 
this paper look visually good, i.e., the least squares procedure 
has converged to an apparently valid solution. However, differ- 
ent initial values might lead to slightly different combinations 
of parameters values beca use, for exampl e, R^.^^ and p, are in- 
trinsically correlated. Malin enet al.l (1201 2l) also noted that small 
changes in the input data can sometimes have a large effect on 
the estimated Plummer parameters. 



^ See |http : //en . wikipedia . org/wiki/Box_plot] 



In models A and B, when observational (i.e., instrumental) 
noise was not included, the parameters derived from N"^™'^ and 
^obs ^gj-g similar to within a few percent. Furthermore, in the 
case of the near distance cloud (beam size equal to 3 pixels) the 
results were not very sensitive to the level of instrumental noise. 
When the resolution of the observations was poor, in particular 
for the far cloud distance of 37 1 pc that corresponded to convolu- 
tion with a beam with FWHM equal to 12 pixels, noise did cause 
significant scatter. In the central density pc the errors could occa- 
sionally increase up to two orders of magnitude. In this case the 
beam is larger than the width of the filament and the observations 
give little constraints on the central profile. Occasionally the re- 
sult is a very narrow and high density peak that, when convolved 
with the large beam, is still consistent with the observations. 

Compared to the model A, the model B was much less af- 
fected by the noise because of its higher column density and 
thus higher signal-to-noise ratios. In model C, with modified 
dust, central density was overestimated by up to ~0.5-l orders 
of magnitude, or even more with the lowest resolution. Similarly, 
the mass estimates could be in error by a factor of several. The 
scatter was much larger than in the models A and B, especially 
for the medium and far cloud distances. For the clouds at the far 
distance (lowest resolution), the total width of the profile is dom- 
inated by the convolution and the values of the derived Plummer 
parameters are more uncertain. This is clearly demonstrated by 
Fig.[TT] The typical errors increase by a factor of several between 
the near and far distances. On the other hand, in the case of the 
model B the errors are only half of those found in the model A. 
For model C the errors are larger. The mass estimates are directly 
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Fig. 9. The FWHM values of the twelve filaments. The frames a, 
b, and c correspond to the models A, B, and C, respectively. The 
squares and the triangles correspond to the near cloud distance 
of 93 pc and the medium cloud distance of 186pc, respectively. 
The solid symbols are values derived from A^t™'= and the open 
symbols values from A^'-"'^ The dashed lines indicate the median 
values for the FWHM values derived from A^°''\ Upper line cor- 
responds to the distance of 186 pc and lower Une to the distance 
of 93pc. 



affected by the wrong assumption of dust opacity but the other 
parameters do not show very significant bias. 

The central question is, are the fitted parameter values, in 
particular p, so reliable that theoretical conclusions can be made 
based on them. In all our models the p parameter values have 
a bipolar behaviour. The lower values {~2-4) are mostly within 
a few percent from those derived from the true column density. 
The largest errors occur when the p value derived from either 
the observed or true column density is near the allowed maxi- 
mum value of 10. In models A and B this occurs only when in- 
strumental noise is included, but in model C also without noise. 
Noise and low resolution clearly increase the probability of get- 
ting a high value for p, which irnplies t hat the parameters are not 
well constrained. Malinen et alJ (12012 ") have shown that the con- 
fidence regions of the Plummer parameters can be wide and the 
parameters (e.g., p and /?flat) strongly correlated. We also com- 
pared the p parameter values derived from our three models for 




4 6 8 10 

Filament number 

Fig. 10. The p values of the twelve filaments using medium 
cloud distance of 186pc (beam size of 6 pixels) and y6 =2.0. 
Squares indicate Model A, triangles Model B, and circles Model 
C. fl) Values are derived from the observed A^'-"'^ without noise. 
b) Values are derived from the observed A^'-''"* with noise, c) 
Values are derived from the true A^^™'^. 



individual filaments (Fig. [TOI i. The effect of noise appears to be 
much stronger than the radiative transfer effects. 

To separate the effects of instrumental noise from those of 
spatial resolution we carried out tests where only the noise lev- 
els were varied. Figure [12] shows the results for model A at the 
near distance. The noise is scaled between 0.5 and 4.0 times the 
values used above. The figure shows the median parameter val- 
ues and their inter-quartile ranges for each filament and each 
noise level, as determined from 50 noise realisations. With the 
default noise, all the parameters shown are still well determined. 
When the noise is increased, the scatter between noise realisa- 
tions increases rapidly. The noise also causes bias in the case of 
some filaments. However, this varies from filament to filament 
and, for example, the median value of Mnne can either increase 
or decrease. The figure also implies that as the cloud distance 
is increased in the simulations, part of the increased scatter in 
the parameter values is caused by noise. At the far distance, the 
filament appears four time smaller than at the near distance and 
the total signal-to-noise ratio (integrated over the source) is cor- 
respondingly smaller. At the relative noise value of 4.0 the value 
of p is well constrained (within 50%) for less than half of the 
selected filaments. In the fits p was always restricted to values 
below 10.0. 

We have concentrated on the differences between the true 
and the observed column densities and their effect on the es- 
timated filament parameters. However, we can note some ad- 
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Fig. 11. Boxplot for the relative errors of Ig pc, p, and Mune for 
the three cloud models and the three cloud distances (beam sizes 
of 3, 6, and 12 pixels corresponding to the near, medium, and 
far distances). The errors are calculated by comparing the values 
derived from A^^*"^ to those obtained from A^^™'' (see Figs.|4]|5] 
and [8ll. For each parameter, the box extends between the first 
and third quartile with a horizontal line indicating the median 
value. The vertical lines (whiskers) extend from the box to the 
last data point within a distance equal to 1.5 times the distance 
between the third and the first quartile. Points outside this range 
are drawn with crosses. The number of outliers falling outside 
the plotted range is indicated with numbers near the upper and 
lower borders. 



ditional factors affecting the interpretation of such parameters. 
Like in our models, the line-of-sight confusion can be a prob- 
lem in real observations. For example, some of the occasional 
large values of the parameter p are suspected to be caused by 
background structures. Those can be confused with the filament 
itself or they may affect the fits by modifying the background 
estimates (see Fig.O. 

The line-of-sight confusion can be even more severe and 
some of the identified filaments can be the result of a chance 
alignment of unrelated (or related but still physically separate) 
density structures. Of the four filaments selected from the x di- 
rection maps, none is seen as a clear isolated filament when ex- 
amined from the other two orthogonal directions. Figure [13] il- 
lustrates this in the case of filament 1. The filament was iden- 
tified from the image towards the direction x. We determined 
the locations of the volume density maxima for a few positions 
along the filament along the line of sight. The images calculated 
for the other two orthogonal viewing directions show that the 
filament is actually a superposition of at least two independent 
density structures. This may be a common scenario also in the 
case of real clouds as suggested, for example, by the identifica- 
tion of several kinematically separate components in some of the 
Taurus filaments (Tafalla 2010, private communication). 

Figure[T4]shows the central densities of filaments as derived 
from the observations. For comparison, the figure includes the 
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Fig. 12. The dependence of the extracted Plummer parameters 
Pc, p, and Miine on the noise level in the case of model A and 
the near cloud distance. The level of instrumental noise is either 
0.5, 1.0, 2.0, or 4.0 times the default value (see Sect. |2]i. The 
plot shows for all the twelve filaments the median values and, as 
vertical bars, the inter-quartile intervals of the estimated param- 
eters that are obtained from 50 noise realisations of the surface 
brightness maps. The coloured lines connect the median values. 
The vertical bars plotted for the individual filaments are shifted 
slightly along the x axis to avoid overlap. 
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Fig. 13. Filament 1 seen from two orthogonal viewing direc- 
tions. The object was identified from observations towards the 
X direction. The figure shows the object from the y and z direc- 
tions (frames a and b, respectively). The circles denote volume 
density maxima for a few lines of sight through the centre of the 
filament as seen towards the x direction. The background image 
is the Ay of the model cloud. 
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actual maximum line-of-sight density that is obtained from the 
model clouds and is averaged along the length of the projected 
filament. For model C the densities are overestimated because 
of the difference between the true and the assumed Ky values. In 
models A and B there is less bias but even in model B the den- 
sities are sometimes overestimated by a factor of a few. This is 
the case for example for the filament number 1 that was already 
seen to consist of separate density structures along the line of 
sight. The result reminds us that continuum observations cannot 
directly measure the volume density. The uncertainty on the pc 
values is several tens of per cent. 

Figure [14] also compares the mass estimates Mw^e with the 
actual cloud mass within a projected distance of 0. 1 pc of the fila- 
ment centre. Apart from model C, the true masses are on average 
slightly higher than the values derived from the fits of the fila- 
ment profiles. To some extent this may reflect the presence of an 
extended cloud component that in the Plummer fits is modelled 
as a separate component and thus is not included in the Miine esti- 
mates. In some cases there is clear anticorrelation between the pc 
and Miine values (e.g., filaments 4 and 8). When the fitted profile 
has a narrow peak, a larger fraction of the mass gets attributed 
to the background component. However, the background is usu- 
ally low (see Fig.O and cannot always explain the factor of ~2 
differences in the models A and B. 

The bias in the mass estimates can be traced back to the as- 
sumptions of the column density calculations, i.e., a constant 
line-of-sight temperature and a dust spectral index yS - 2.0 
(instead of the slightly larger actual value of the dust model). 
Although the bias in the column density values should increase 
towards the centre of the filament, this was not reflected in the 
parameters pc or /? as a systematic error In model C the inclusion 
of dust with higher sub-millimetre opacity causes bias mainly in 
the mass estimates but also the central density and the parame- 
ter p should be affected. The fact that the column densities were 
overestimated more towards the centre of the filaments should 
lead to overestimation of the parameter pc- This is also seen in 
the results, the central densities being overestimated by a factor 
of 2-3. This is consistent with the modified dust being respon- 
sible for most of the signal along the central line of sight. The 
relative abundance of the dust with higher Ky was increased as 
a function of density and this could have been expected to lead 
to higher p values. In practice, no clear bias is seen. Together 
with the high column density, the high dust emissivity increased 
the signal-to-noise ratio in model C far above that of models A 
and B. Nevertheless, in model C also the scatter of the parameter 
values was the largest. 



6. Conclusions 

We have examined the reliability of the parameters that can be 
derived for interstellar cloud filaments using dust emission data. 
The study is based on synthetic observations calculated with 
MHD simulations of filamentary clouds and with radiative trans- 
fer modelling of the sub-millimetre dust emission. The observa- 
tional parameters correspond to typical Herschel observations. 
We have examined the estimated FWHM widths of the filaments 
and the parameters of the Plummer functions used to describe 
the filament profiles. The results have led to the following con- 
clusions: 

- For clouds at a distance of -lOOpc, the profile parameters 
obtained from the observations and from the true column 
density are usually similar to within a few percent. However, 
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Fig. 14. Comparison of the pc and filament mass estimates from 
the Plummer fits (open symbols) and the actual values of maxi- 
mum line-of-sight densities and masses in the cloud model (solid 
squares). The error bars on the true densities correspond to the 
1 cr dispersion of the peak density values along the filament. In 
the mass estimates the solid squares correspond to the total mass 
within a projected distance of 0. 1 pc of the filament. Model A: 
frames a and d, model B: frames b and e, model C: frames c and 
/■ 



differences up to tens of percent are observed in the case of 
a few filaments. 

- The results are affected by the instrumental noise. In the 
higher density models the signal-to-noise ratios were higher 
and the errors were correspondingly smaller If the instru- 
mental noise is increased from the assumed default values, 
the errors increase rapidly. When the noise is increased by 
a factor of four, the filament parameters are no longer well 
constrained even for nearby clouds. 

- The uncertainties caused by the background confusion and 
the typical instrumental noise outweigh the radiative transfer 
effects. 

- The uncertainties increase rapidly when the resolution of the 
observations is decreased. At a distance of ~400pc the fila- 
ments are barely resolved, the typical errors are larger by a 
factor of a few, and there are more cases where the profile fit 
leads to completely wrong parameter values. 

- We examined a model where the dust opacity is increased as 
a function of density. In the analysis the mass estimates are 
naturally affected by wrong assumptions of the dust opacity. 
The errors are larger also for the other parameters but they 
are not strongly biased. 

- The derived filament properties are affected by the line-of- 
sight confusion. In our study some of the filaments that were 
identified from the column density maps are not continuous 
structures of the three dimensional density field. 
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